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Bifurcations without Parameters: 
Some ODE and PDE Examples 



Bernold Fiedler* Stefan Liebscher^ 

Abstract 

Standard bifurcation theory is concerned with families of vector fields 

x = f(x, A) 

x G R™, involving one or several constant real parameters A. Viewed as a 
differential equation for the pair (x,X), we observe a foliation of the total 
phase space by constant A. Frequently, the presence of a trivial stationary 
solution x — is also imposed: 

= /(0,A). 

Bifurcation without parameters, in contrast, discards the foliation by a 
constant parameter A. Instead, we consider systems 

* = /(a>> V) 

y = g{x,y) 

Standard bifurcation theory then corresponds to the special case y = A, g = 0. 
To preserve only the trivial solution x = 0, instead, we only require 

o = mv)=g(o,y) 

for all y. A rich dynamic phenomenology arises, when normal hyperbolicity 
of the trivial stationary manifold x — fails, due to zero or purely imaginary 
eigenvalues of the Jacobian fx{0,y)- 

Specifically, we address the cases of failure of normal hyperbolicity due 
to a simple eigenvalue zero, a simple purely imaginary pair (Hopf bifurcation 
without parameters), a double eigenvalue zero (Takens-Bogdanov bifurcation 
without parameters) , and due to a double eigenvalue zero with additional time 
reversal symmetries. 

We illustrate consequences of our results with examples from ordinary and 
partial differential equations arising in systems of coupled oscillators, in the 
analysis and numerics of hyperbolic conservation and balance laws, and in the 
fluid dynamics of plane Kolmogorov flows. 

The results are joint work with Andrei Afendikov, James C. Alexander, 
and Stefan Liebscher. For references see 
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1. Applied motivation 

In this article we sketch and illustrate some elements of the nonlinear dynamics 
near equilibrium manifolds. Denoting the equilibrium manifold by x = 0, in local 
coordinates (x, y) £ IR n X IR fe , we consider systems 

i = f(x,y) 

y = 9\ x ,y) 

and assume 

/(0,y) = g(0,y) = 0, (1.2) 

for all y. For simplicity we will only address the cases k = 1 of lines of equilibria, 
and k = 2 of equilibrium planes. Sufficient smoothness of /, g is assumed. The 
occurrence of equilibrium manifolds is infinitely degenerate, of course, in the space 
of all vector fields (/, g) - quite like many mathematical structures are: equivariancc 
under symmetry groups, conservation laws, integrability, symplectic structures, and 
many others. The special case 

g ee (1.3) 

in fact amounts to standard bifurcation theory, in the presence of a trivial solution 
x — 0; see for example [5]. Note that condition l|1.3fl . which turns the fc-dimensional 
variable y into a preserved constant parameter, is infinitely degenerate even in 
our present setting (|1.2|) of equilibrium manifolds. Due to the analogies of our 
results and methods with bifurcation theory, we call our emerging theory bifurcation 
without parameters. This terminology emphasizes the intricate dynamics which 
arises when normal hyperbolicity of the equilibrium manifold fails; see the sections 
below. 

To motivate assumption (|1.2(l . we present several examples. First, consider 
an "octahedral" graph T of 2(m + 1) vertices {±1, . . . , ±(m + 1)}. The graph T 
results from the complete graph by eliminating the "diagonal" edges, which join the 
antipodal vertices ±jf, for j = 1, . . . , m + 1. For m — 1 we obtain the square, for 
m = 2 the octahedron, and so on. Consider the system 

u>j = fj(uj, Ufe ) (1.4) 

of oscillators Uj £ \R n on T, additively coupled along the edges by fj. We assume 
an antipodal oddness symmetry of the individual oscillator dynamics 



f-j{-Uj,0) = -fj(u j} 0). 



(1.5) 
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As a consequence, the antipode space 

£ := {u=(uj) je r; U-j = -Uj} (1.6) 

is invariant under the flow 1)1. 4[l . Moreover, the flow on £ completely decouples into 
a direct product flow of the m + 1 diagonally antipodal, decoupled pairs 

u±j = f±j(u ±j ,0). (1.7) 

For the square case m = 1, this decoupling phenomenon was first observed in 
For more examples see also 

An m-plane of equilibria arises from periodic solutions of the decoupled system 
(|1.7fl . Assume (II .7|) possesses time periodic orbits Uj(t + ipj) of equal period Tj — 
2tt, for j = l,...,m + 1. Choose arbitrary phases ^ € 5 1 and let U-j(t) :— 
—Uj(t), tp-j — ifij. Then 

:= ( % -(t + ^-)W G £ (1.8) 

is a 27r-periodic solution of 11.4|l . 1)1.71) . for arbitrary phases ip £ T m+1 . Eliminating 
one phase angle ip m +i by passing to an associated Poincare map, an m-dimensional 
manifold of fixed points arises, parametrized by the remaining m phase angles. 
Assuming, in addition to the diagonal oddness symmetry l|1.5ll . equivariance of (|1.4|l 
with respect to an 5 1 -action, the Poincare map can in fact be obtained as a time-27r 
map of an autonomous flow within the Poincare section. In suitable notation, y = 
{<£>!, . . . , v?m}, the fixed point manifold then becomes an m-dimensional manifold of 
equilibria, as presented in l|l.l|l . 1)1. 2|l above. For more detailed discussions of this 
example in the context of bifurcations without parameters see EH ■ 

As a second example of equilibrium manifolds we consider viscous profiles 
u = — s t)l £ ) °f systems of nonlinear hyperbolic conservation laws and stiff 
balance laws 

d t u + d^F{u) +e- 1 G(u) = ed^u. (1.9) 
Viscous profiles then have to satisfy an e-independent ODE system 

u = (F'(u)- s-id)u + G(u). (1.10) 

Standard conservation laws, for example, require G = 0. The presence of m conser- 
vation laws corresponds to nonlinearities G with range in a manifold of codimension 
m in u-space. Typically, then, G(u) = describes an equilibrium manifold of di- 
mension m of pairs (u, ii) — (u, 0), in the phase space of 1)1.10)1 . For an analysis of 
this example in the context of bifurcations without parameters see |H1 116). For an- 
other example, which relates binary oscillations in central difference discretizations 
of hyperbolic balance laws with diagonal uncoupling of coupled oscillators, see jll| . 

We conclude our introductory excursion with a brief summary of some further 
examples. In [7], lines of equilibria have been observed for the dynamics of models 
of competing populations. This included a first partial analysis of failure of normal 
hyperbolicity. 

A topologically very interesting example in compact three-dimensional mani- 
folds involves contact structures r/(^) (i.e., nonintegrablc plane fields and gradient 
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vector fields £ 



VV(£) € ??(£)• See jS] for an in-depth analysis. Examples include 



mechanical systems with nonholonomic constraints. Notably, level surfaces of reg- 
ular values of the potential V consist of tori. Under a nondegeneracy assumption, 
equilibria form embedded circles, that is, possibly linked and nontrivial knots. 

For a detailed study of plane Kolmogorov fluid flows in the presence of a 
line of equilibria with a degeneracy of Takens-Bogdanov type and an additional 
reversibility symmetry, see pQ. 

As a caveat we repeat that lines of equilibria, which are transverse to level 
surfaces of preserved quantities A do not provide bifurcations, without parameters. 
Rather, y = for y := A exhibits this problem as belonging to standard bifurcation 
theory; sec HI. 311 . 

2. Sample vector fields and resulting flows 

In this section we collect relevant example vector fields (11.1(1 . I|1.2(l with lines 
and planes of equilibria x = 0; see ^3 for further details. We illustrate and 
comment the resulting flows. 

Normally hyperbolic equilibrium manifolds admit a transverse C°-foliation 
with hyperbolic linear flows in the leaves. See for example [TS], [201 an d the ample 
discussion in 0|. As a first example, we therefore consider 



Note the loss of normal hyperbolicity at x = y = 0, due to a nontrivial transverse 
eigenvalue zero of the linearization. Clearly dx/dy = y, and the resulting flow lines 
are parabolas; see Figure^ For comparison with standard bifurcation theory, where 
y = A, we draw the j/-axis of equilibria horizontally. 
As a second example, consider 



Again, a transverse zero eigenvalue occurs - this time with an additional reflection 
symmetry y » —y. Dividing by the Euler multiplier x, the reflection becomes a 
time reversibility. See the left parts of Figure [21 for the resulting flows. Note the 
resulting integrable, harmonic oscillator case which originates from the elliptic sign 
y = -x 2 . 

As a third example, we consider x — (xi, X2) € IR 2 , y € IR with a line x = 
of equilibria and a purely imaginary nonzero eigenvalue iui at x — 0. Normal- 
form theory, for example as in |21| . then generates an additional ^-symmetry by 
the action of exp(iwt) in the x-eigenspace. This equivariance can be achieved, 
successively, up to Taylor expansions of any finite order. In polar coordinates (r, (p) 
for x, an example of leading order terms is given by 



x = xy, 
y = x. 



(2.1) 



x = xy, 
y = ±x 2 . 



(2.2) 



r = ry, 




(2.3) 
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Since the first two equations in (|2.3[) coincide with l|2.2[) . the dynamics is then 
obtained by simply rotating the left parts of Figure |2 around the y-axis at speed 
w. The right parts of Figure [5] provide three-dimensional views of the effects of 
higher-order terms which do not respect the ^-symmetry of the normal forms. In 
the elliptic case (b) , all nonstationary orbits are heteroclinic from the unstable foci, 
at y > 0, to stable foci at y < 0. The two-dimensional respective strong stable and 
unstable manifolds will split, generically, to permit transverse intersections. 

Our fourth example addresses Takens-Bogdanov bifurcations without param- 
eters. In suitable rescaled form it reads 

y + yy = e((A - y)y + by 2 ) (2.4) 

with fixed parameters b, A and e. The y-axis as equilibrium line is complemented 
by the two transverse directions x — (y, y). Note the algebraically triple zero 
eigenvalue, double in the transverse ^-directions, for A = y = 0. Two examples of 
the resulting dynamics for small positive e are summarized in Figure |3J 

The coordinates r and H in Figure|3]are adapted to the completely integrable 
case e = 0. Indeed, obvious first integrals are then given by O = y + \y 2 and 
H = \y 2 — yy — -|y 3 . Coordinates are r = logG and H = 0~ 3 / 2 H, not drawn 
to scale. Parameters are e, A > and, for the hyperbolic case, —17/12 < b < — 1. 
For the elliptic case we consider b > — 1. The equilibrium y-axis, a cusp in (0,i?) 
coordinates, transforms to the top (saddles) and bottom (foci) horizontal boundaries 
H = ±|\/2, with y = relegated to r = — oo. Since r and H are constants of the 
flow, for e = 0, they represent slow drifts on the unperturbed periodic motion, for 
small e > and | H |< |\/2- The top value H = +|V2 also represents homoclinics 
to the saddles, for e — 0. 
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Figure 2: Lines of equilibria (y-axis) with imaginary eigenvalues: 
Hopf bifurcation without parameters. Case (a) hyperbolic; case 
(b) elliptic. Red: strong unstable manifolds; green: stable mani- 
folds. 




Figure 3: Takens-Bogdanov bifurcations without parameters. Case (a) hyperbolic; case (b) elliptic. 
Top: stable and unstable manifolds; bottom: invariant sets. For coordinates and fixed parameters 
see text. 
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Along the focus line H = —\\F^ we observe Hopf bifurcations without pa- 
rameters, corresponding to y = X > 0. The value of b distinguishes elliptic and 
hyperbolic cases. In addition, lines of saddle-focus heteroclinic orbits and isolated 
saddle-saddle heteroclinics are generated, for e > 0, by breaking the homoclinic 
sheets of the integrable case. Note in particular the infinite swarm of saddle-saddle 
heteroclinics, in the hyperbolic case. 

As a final, fifth example we consider a reversible Takens-Bogdanov bifurcation 
without parameters: 

y+(l-3y 2 )y = ayij + by 2 . (2.5) 

Here we fix a, b to be small. Again x — (y, y) denotes the directions transverse to 
the equilibrium y-axis. Time reversibility generates solutions —y(—t) from solutions 
y(t). For two examples of the resulting dynamics see Figure Coordinates are the 
obvious first integrals, for e — 0, given by = y + y — y 3 and H = —yy + -|y 2 + 
|y 4 — \y 2 ■ Note the two Takens-Bogdanov cusps, separated by a Hopf point along 
the lower arc of the equilibrium "triangle". Compare Figures|5J|21 The elliptic Hopf 
point (b) arises for a ■ (a — b) > 0, whereas a - (a — b) < in the hyperbolic case (a). 
Also note the associated finite and infinite swarms of saddle-saddle heteroclinics, 
respectively. 

3. Methods 

Pictures are not proofs. What has been proved, then, and how? We use 
ingredients involving algebra, analysis, and numerical analysis, as we outline in this 
section. For further details see |16l 1101 151 H). 

In a first algebraic step, we derive normal forms for vector fields with lines 
or planes of equilibria, assuming spectral degeneracies of the linearization A in 
transverse directions. The spectral assumptions on A in fact coincide with those 
established for parameter-dependent matrix families in standard bifurcation theory. 
This is reflected in the naming of the five examples of section [21 

There are more or less standard procedures to derive normal forms of vec- 
tor fields. By suitable polynomial diffeomorphisms, certain Taylor coefficients of 
the vector field are eliminated, successively, for higher and higher order. See for 
example |21| for a systematic choice of normal forms, particularly apt for introduc- 
ing equivariance of the nonlinear normal- form terms under the action of ex~p(A T t). 
Normal forms are, however, nonunique in general and other choices are possible. 

In the present cases, we adapt the normal-form procedure to preserve the lo- 
cally flattened equilibrium manifolds. Although the approach in |21| can be modified 
to accommodate that requirement, it did not provide vector fields convenient for 
subsequent analysis of the flow. A systematic approach to this combined problem is 
not known, at present. All examples l|2.1|l - H2.5(l represent truncated normal forms. 
For the derivation of specific normal forms, for example of 12.4J1 . to any order, and 
of example Ij2.5|l . to third order, see 0,0, respectively. 

Subsequent analysis of the normal- form vector fields is based on scalings, alias 
blow-up constructions. This is the origin, for example, of the small scaling pa- 
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Figure 4: Reversible Takens-Bogdanov bifurcations without pa- 
rameters. Case (a) hyperbolic; case (b) elliptic. For coordinates 
and fixed parameters see text. 
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rameter e in the Takens-Bogdanov example (|2.4I) . In passing, we note a curious 
coincidence of two view points for (|2.4() . concerning the roles of the equilibrium 
coordinate y S IR and the "fixed" real parameter A G IR. First, we may consider A 
as a parameter, with a line of equilibria y associated to each fixed A. Then (|2.4|l 
describes the collision of a transverse zero eigenvalue at y = 0, as in l|2.1|L with 
imaginary Hopf eigenvalues at y = A > 0, as in (|2.3|l . as A decreases through zero. 
Alternatively, we may consider normal forms for a plane y = (2/1,2/2) of equilibria 
with a transverse double zero eigenvalue, at y = 0. It turns out that the two cases 
coincide, after a scaling blow-up, up to second order in e, via the correspondence 
V =2/i>A = 2/2. 

The core of any successful flow analysis in bifurcation theory is an integrable 
vector field; see again sectional The issues of nonintegrable perturbations, by small 
e > 0, and of omitted higher order terms, not in normal form, both ensue. Since 
the underlying integrable dynamics is periodic or homoclinic, in examples (|2.3|l - 
<|2.5[) . averaging procedures apply. Indeed, e > then introduces a periodically 
forced, slow flow on first integrals, like (0, H), characteristic of e = 0. We therefore 
derive an appropriate, but autonomous Poincare flows, on (0,i?), such that the 
associated true Poincare map can be viewed as a time discretization of first order 
and step size e. In the unperturbed periodic region, this amounts to averaging, 
while the Poincare flow indicates Melnikov functions at homoclinic or heteroclinic 
boundaries. The exponential averaging results by Neishtadt l^j, for example, then 
imply that the separatrix splittings in the elliptic Hopf case (b) , indicated in Figure 
121 are exponentially small in the radius r of the split sphere, for analytic vector 
fields. See also |T2|. 

Lower bounds of separatrix splittings have not been established, in our set- 
tings. This problem is related to the very demanding Lazutkin program of asymp- 
totic expansions for exponentially small separatrix splittings. For recent progress, 
including the case of Takens-Bogdanov bifurcations for analytic maps, see ^B] an d 
the references there. In absence of rigorous lower bounds, our figures indicate only 
simplest possible splitting scenarios. 

While the splitting near elliptic Hopf points are exponentially small, the dis- 
cretization of the Poincare flow also exhibits splittings of the unperturbed saddle 
homoclinic families, which are of first order in the perturbation parameter e, in ex- 
ample (I2.4() . or in the small parameters a, 6, in example (|2.5|l . Explicit expressions 
have been derived for the Melnikov functions associated to these homoclinic split- 
tings, in terms of elliptic function in case |23, and even of elementary functions in 
case Ij2.5|l . Simplicity and uniqueness of zeros of the Melnikov functions, however, 
has only been confirmed numerically. While this does not, strictly speaking, match 
an analytic proof, it still at least supports the validity of the scenarios summarized 
in Figures |3 and 0] 

4. Interpretation and perspective 

We indicate some consequences of the above results for the examples of coupled 
oscillators, viscous shock profiles, and Kolmogorov flows indicated in section ^ 
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We conclude with a few remarks on the future perspective of bifurcations without 
parameters. 

We first return to the example (|1.4[l - (|1.8f) of a coupled oscillator square, m = 1. 
Equilibria y of the Poincare flow then indicate decoupled antipodal periodic pairs, 
say with phase difference y + c. The case of a transverse zero eigenvalue, l|2.1[) and 
Figure n then indicates a 50% chance of recovery of decoupling with a stable phase 
difference y < 0, locally, even when the stability threshold y = has been exceeded. 
The hyperbolic Hopf case (|2.2|) . Figure El (a), illustrates immediate oscillatory loss 
of decoupling stability by transverse imaginary eigenvalues. A 100% recovery of 
decoupling stability, in contrast, occurs at elliptic Hopf points; see H2.2(l . Figure El 
(b). The exponentially small Neishtadt splitting of separatrices indicates a very 
delicate variability in the asymptotic phase relations of this recovery, for t — > ±oo. 
See UDJ. The Takens-Bogdanov cases (|2.4|l . Figures |3] (a), (b) can then be viewed 
as consequences of a mutual interaction, of a transverse zero eigenvalue with cither 
Hopf case, for recovery of stable decoupling. 

In the example (|1.9|) . (|1.10|) of stiff balance laws, elliptic Hopf bifurcation 
without parameters as in (|2.4|) , Figure [3| (b) , indicates oscillatory shock profiles 
u(r), t = (£ — st)/e. Such profiles in fact contradict the Lax condition, being over- 
compressive, and violate standard monotonicity criteria. For small viscosities e > 0, 
weak viscous shocks in fact turn out unstable, in any exponentially weighted norm, 
unless they travel at speeds s exceeding all characteristic speeds. The oscillatory 
profiles can be generated, in fact, by the interaction of inherently non-oscillatory 
gradient flux functions F(u) with inherently non-oscillatory gradient-like kinetics 
G(u) in systems of dim > 3. See [Ifi] . 

The problem of plane stationary Kolmogorov flows asks for stationary solutions 
of the incompressible Navier-Stokes equations in a strip domain (£, rj) £ IR x [0, 2ir], 
under periodic boundary conditions in 77; see |17j . An 77-periodic external force 
(F(r]),0), is imposed, acting in the unbounded ^-direction. Kolmogorov chose 
F(rj) = sin 77. The Kirchgassner reduction 14 captures all bounded solutions which 
are nearly homogeneous in £, in a center manifold spirit which lets us interpret £ as 
"time". The resulting ordinary differential equations in IR 6 reduce to IR 3 , by fixing 
the values of three first integrals. A line of ^-homogeneous equilibria appears, in 
fact, and Kolmogorov's choice corresponds to example (|2.5|) with a = b = 0, to lead- 
ing orders. In particular note the double reversibility, then, under y(t) t— > ±y(—t) 
which is generated by 



As observed by Kolmogorov, an abundance of spatially periodic profiles results. 
The sample choice F{rj) — sin 77 + csin2?7, in contrast, which breaks the first of the 
symmetries in (I4.1|l , leads to (|2.5|) with b = < a, alias an elliptic reversible Takens- 
Bogdanov point without parameters; see Figure^] (b). In particular, the set of near- 
homogeneous bounded velocity profiles of the incompressible, stationary Navier- 
Stokes system is then characterized by an abundance of oscillatory heteroclinic 
wave fronts, which decay to different asymptotically homogeneous (^-profiles, for 
£ — * ±00. The PDE stability of these heteroclinic profiles is of course wide open. 



F( v ) 



F(rj + ir), 
F(-rj). 



and 



(4.1) 
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As for perspectives of our approach, we believe to have examples at hand, from 
sufficiently diverse origin, to justify further development of a theory of bifurcations 
without parameters. In fact, transverse spectra {0,±iu} and {±iwi,±iw 2 } still 
await investigation before we can claim any insight into nonhyperbolicity of even 
the simple case of an equilibrium plane. This assumes the absence of further struc- 
tural ingredients like symplecticity, contact structures, symmetries, and the like. 
Certainly our example collecting activities are far from complete, at this stage. 

In addition, we have not addressed the issue of perturbations, so far, which 
could destroy the equilibrium manifolds by small drift terms. Examples arise, for 
example, when slightly detuning the basic frequencies 2ir/Tj of our uncoupled oscil- 
lators or, much more generally, in the context of multiple scale singular perturbation 
problems. Feedback and input from our readers will certainly be most appreciated! 
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